Geographic variability of hybridization between red-breasted and red-naped sapsuckers.

Hybrid zones reveal the strength of reproductive isolation between populations undergoing speciation and are a key tool in evolutionary biology research. Multiple replicate transects across the same hybrid zone offer insight into the dynamics of hybridization in different environments, clarifying the role of extrinsic forces on the speciation process. Red-breasted and red-naped sapsuckers (Sphyrapicus ruber and Sphyrapicus nuchalis) have a long zone of contact over approximately 1,600 km from central British Columbia, Canada to central California, USA. We used Genotyping-by-Sequencing data from three independent sapsucker hybrid zone transects to compare hybridization dynamics between these species under variable geoclimatic conditions. We generated geographic clines of the genomic data to compare hybrid zone widths and used random forests models and linear regression to assess the relationship between climate and sapsucker ancestry along each transect. Our results show variation in the directionality of backcrossing, often indicative of moving hybrid zones. We note variable cline widths among transects, indicating differences in selection maintaining hybrid zone dynamics. Furthermore, random forests models identified different variables in close association with sapsucker ancestry across each transect. These results indicate a lack of repeatability across replicate transects and a strong influence of the local environment on hybrid zone dynamics.

Estimating hybridization in the wild using citizen science data: A path forward.

Genomic evidence of introgression in natural populations has reinvigorated the study of hybridization in recent years. Still, it is largely unknown how frequently individual organisms mate across species lines. Recently, Justyn et al. suggested that eBird, one of the world's largest citizen science databases, may supply adequate data for estimating hybridization rates. Here, we compare Justyn et al.'s estimates-and their conclusions that hybridization is rare-with estimates from museum and molecular data. We also estimate hybridization using eBird observations from areas and times when hybridization is possible, namely, in contact zones during the breeding season. These estimates are all considerably higher than those reported in Justyn et al., emphasizing that inferences from multiple datasets can differ radically. Finally, we demonstrate an approach for predicting the location of hybrid zones using eBird data, which can be done with high confidence and with unprecedented resolution. We show that citizen science data, far from settling the question of how frequently bird species hybridize, instead offer a promising step toward more focused study of hybrid zones.

Temperature heterogeneity correlates with intraspecific variation in physiological flexibility in a small endotherm.

Phenotypic flexibility allows individuals to reversibly modify trait values and theory predicts an individual's relative degree of flexibility positively correlates with the environmental heterogeneity it experiences. We test this prediction by integrating surveys of population genetic and physiological variation with thermal acclimation experiments and indices of environmental heterogeneity in the Dark-eyed Junco (Junco hyemalis) and its congeners. We combine field measures of thermogenic capacity for 335 individuals, 22,006 single nucleotide polymorphisms genotyped in 181 individuals, and laboratory acclimations replicated on five populations. We show that Junco populations: (1) differ in their thermogenic responses to temperature variation in the field; (2) harbor allelic variation that also correlates with temperature heterogeneity; and (3) exhibit intra-specific variation in thermogenic flexibility in the laboratory that correlates with the heterogeneity of their native thermal environment. These results provide comprehensive support that phenotypic flexibility corresponds with environmental heterogeneity and highlight its importance for coping with environmental change.

Phenotypic and genetic introgression across a moving woodpecker hybrid zone.

In hybrid zones in which two divergent taxa come into secondary contact and interbreed, selection can maintain phenotypic diversity despite widespread genetic introgression. Red-breasted (Sphyrapicus ruber) and red-naped (S. nuchalis) sapsuckers meet and hybridize along a narrow contact zone that stretches from northern California to southern British Columbia. We found strong evidence for changes in the structure of this hybrid zone across time, with significant temporal shifts in allele frequencies and in the proportions of parental phenotypes across the landscape. In addition to these shifts, we found that differences in plumage predict genetic differences (R2  = 0.80), suggesting that plumage is a useful proxy for assessing ancestry. We also found a significant bimodal distribution of hybrids across the contact zone, suggesting that premating barriers may be driving reproductive isolation, perhaps as a result of assortative mating based on plumage differences. However, despite evidence of selection and strong patterns of population structure between parental samples, we found only weak patterns of genetic divergence. Using museum specimens and genomic data, this study of sapsuckers provides insight into the ways in which phenotypic and genetic structure have changed over a 40-year period, as well as insight into the mechanisms that may contribute to the maintenance of the hybrid zone over time.

Origin and cross-century dynamics of an avian hybrid zone.

BACKGROUND: Characterizations of the dynamics of hybrid zones in space and time can give insights about traits and processes important in population divergence and speciation. We characterized a hybrid zone between tanagers in the genus Ramphocelus (Aves, Thraupidae) located in southwestern Colombia. We evaluated whether this hybrid zone originated as a result of secondary contact or of primary differentiation, and described its dynamics across time using spatial analyses of molecular, morphological, and coloration data in combination with paleodistribution modeling. RESULTS: Models of potential historical distributions based on climatic data and genetic signatures of demographic expansion suggested that the hybrid zone likely originated following secondary contact between populations that expanded their ranges out of isolated areas in the Quaternary. Concordant patterns of variation in phenotypic characters across the hybrid zone and its narrow extent are suggestive of a tension zone, maintained by a balance between dispersal and selection against hybrids. Estimates of phenotypic cline parameters obtained using specimens collected over nearly a century revealed that, in recent decades, the zone appears to have moved to the east and to higher elevations, and may have become narrower. Genetic variation was not clearly structured along the hybrid zone, but comparisons between historical and contemporary specimens suggested that temporal changes in its genetic makeup may also have occurred. CONCLUSIONS: Our data suggest that the hybrid zone likey resulted from secondary contact between populations. The observed changes in the hybrid zone may be a result of sexual selection, asymmetric gene flow, or environmental change.

Changing climate mediates sapsucker (Aves: Sphyrapicus) hybrid zone movement.

Hybrid zones, where two divergent taxa meet and interbreed, offer unique opportunities to investigate how climate contributes to reproductive isolation between closely related taxa and how these taxa may respond to climatic changes. Red-naped (Sphyrapicus nuchalis) and Red-breasted (Sphyrapicus ruber) sapsuckers (Aves: Picidae) hybridize along a narrow contact zone that stretches from northern California to British Columbia. The hybrid zone between these species has been studied extensively for more than 100 years and represents an excellent system for investigations of the evolution of reproductive isolation. Shifts in the proportions of phenotypes at hybrid localities since 1910 that were inferred using specimens from museum collections were confirmed using species distribution models. We predicted the historical, current, and future distributions of parental and hybrid sapsuckers using Random Forests models to quantify how climate change is affecting hybrid zone movement in the Pacific Northwest. We found observed distribution shifts of parental sapsuckers were likely the result of climate change over the past 100 years, with these shifts predicted to continue for both sapsuckers over the next 80 years. We found Red-breasted Sapsuckers are predicted to continue to expand, while Red-naped Sapsuckers are predicted to contract substantially under future climate scenarios. As a result of the predicted changes, the amount of overlap in the distribution of these sapsuckers may decrease. Using hybrid phenotypes, we found the climate niche occupied by the hybrid zone is predicted to disappear under future conditions. The disappearance of this climate niche where the two parental species come into contact and hybridize may lead to a substantial reduction in genetic introgression. Understanding the impacts of global climate change on hybrid zones may help us to better understand how speciation has been shaped by climate in the past, as well as how evolution may respond to climate change in the future.

Spatially variable coevolution between a haemosporidian parasite and the MHC of a widely distributed passerine.

The environment shapes host-parasite interactions, but how environmental variation affects the diversity and composition of parasite-defense genes of hosts is unresolved. In vertebrates, the highly variable major histocompatibility complex (MHC) gene family plays an essential role in the adaptive immune system by recognizing pathogen infection and initiating the cellular immune response. Investigating MHC-parasite associations across heterogeneous landscapes may elucidate the role of spatially fluctuating selection in the maintenance of high levels of genetic variation at the MHC. We studied patterns of association between an avian haemosporidian blood parasite and the MHC of rufous-collared sparrows (Zonotrichia capensis) that inhabit environments with widely varying haemosporidian infection prevalence in the Peruvian Andes. MHC diversity peaked in populations with high infection prevalence, although intra-individual MHC diversity was not associated with infection status. MHC nucleotide and protein sequences associated with infection absence tended to be rare, consistent with negative frequency-dependent selection. We found an MHC variant associated with a ∽26% decrease in infection probability at middle elevations (1501-3100 m) where prevalence was highest. Several other variants were associated with a significant increase in infection probability in low haemosporidian prevalence environments, which can be interpreted as susceptibility or quantitative resistance. Our study highlights important challenges in understanding MHC evolution in natural systems, but may point to a role of negative frequency-dependent selection and fluctuating spatial selection in the evolution of Z. capensisMHC.

Variation in positively selected major histocompatibility complex class I loci in rufous-collared sparrows (Zonotrichia capensis).

The major histocompatibility complex (MHC) is a highly variable family of genes involved in parasite recognition and the initiation of adaptive immune system responses. Variation in MHC loci is maintained primarily through parasite-mediated selection or disassortative mate choice. To characterize MHC diversity of rufous-collared sparrows (Zonotrichia capensis), an abundant South American passerine, we examined allelic and nucleotide variation in MHC class I exon 3 using pyrosequencing. Exon 3 comprises a substantial portion of the peptide-binding region (PBR) of class I MHC and thus plays an important role in intracellular pathogen defense. We identified 98 putatively functional alleles that produce 56 unique protein sequences across at least 6 paralogous loci. Allelic diversity per individual and exon-wide nucleotide diversity were relatively low; however, we found specific amino acid positions with high nucleotide diversity and signatures of positive selection (elevated d N /d S ) that may correspond to the PBR. Based on the variation in physicochemical properties of amino acids at these "positively selected sites," we identified ten functional MHC supertypes. Spatial variation in nucleotide diversity and the number of MHC alleles, proteins, and supertypes per individual suggests that environmental heterogeneity may affect patterns of MHC diversity. Furthermore, populations with high MHC diversity have higher prevalence of avian malaria, consistent with parasite-mediated selection on MHC. Together, these results provide a framework for subsequent investigations of selective agents acting on MHC in Z. capensis.

Integrating evolutionary and functional tests of adaptive hypotheses: a case study of altitudinal differentiation in hemoglobin function in an Andean Sparrow, Zonotrichia capensis.

In air-breathing vertebrates, the physiologically optimal blood-O2 affinity is jointly determined by the prevailing partial pressure of atmospheric O2, the efficacy of pulmonary O2 transfer, and internal metabolic demands. Consequently, genetic variation in the oxygenation properties of hemoglobin (Hb) may be subject to spatially varying selection in species with broad elevational distributions. Here we report the results of a combined functional and evolutionary analysis of Hb polymorphism in the rufous-collared sparrow (Zonotrichia capensis), a species that is continuously distributed across a steep elevational gradient on the Pacific slope of the Peruvian Andes. We integrated a population genomic analysis that included all postnatally expressed Hb genes with functional studies of naturally occurring Hb variants, as well as recombinant Hb (rHb) mutants that were engineered through site-directed mutagenesis. We identified three clinally varying amino acid polymorphisms: Two in the α(A)-globin gene, which encodes the α-chain subunits of the major HbA isoform, and one in the α(D)-globin gene, which encodes the α-chain subunits of the minor HbD isoform. We then constructed and experimentally tested single- and double-mutant rHbs representing each of the alternative α(A)-globin genotypes that predominate at different elevations. Although the locus-specific patterns of altitudinal differentiation suggested a history of spatially varying selection acting on Hb polymorphism, the experimental tests demonstrated that the observed amino acid mutations have no discernible effect on respiratory properties of the HbA or HbD isoforms. These results highlight the importance of experimentally validating the hypothesized effects of genetic changes in protein function to avoid the pitfalls of adaptive storytelling.

The influence of sampling design on species tree inference: a new relationship for the New World chickadees (Aves: Poecile).

In this study, we explore the long-standing issue of how many loci are needed to infer accurate phylogenetic relationships, and whether loci with particular attributes (e.g., parsimony informativeness, variability, gene tree resolution) outperform others. To do so, we use an empirical data set consisting of the seven species of chickadees (Aves: Paridae), an analytically tractable, recently diverged group, and well-studied ecologically but lacking a nuclear phylogeny. We estimate relationships using 40 nuclear loci and mitochondrial DNA using four coalescent-based species tree inference methods (BEST, *BEAST, STEM, STELLS). Collectively, our analyses contrast with previous studies and support a sister relationship between the Black-capped and Carolina Chickadee, two superficially similar species that hybridize along a long zone of contact. Gene flow is a potential source of conflict between nuclear and mitochondrial gene trees, yet we find a significant, albeit low, signal of gene flow. Our results suggest that relatively few loci with high information content may be sufficient for estimating an accurate species tree, but that substantially more loci are necessary for accurate parameter estimation. We provide an empirical reference point for researchers designing sampling protocols with the purpose of inferring phylogenies and population parameters of closely related taxa.

Species trees for the tree swallows (Genus Tachycineta): an alternative phylogenetic hypothesis to the mitochondrial gene tree.

The New World swallow genus Tachycineta comprises nine species that collectively have a wide geographic distribution and remarkable variation both within- and among-species in ecologically important traits. Existing phylogenetic hypotheses for Tachycineta are based on mitochondrial DNA sequences, thus they provide estimates of a single gene tree. In this study we sequenced multiple individuals from each species at 16 nuclear intron loci. We used gene concatenated approaches (Bayesian and maximum likelihood) as well as coalescent-based species tree inference to reconstruct phylogenetic relationships of the genus. We examined the concordance and conflict between the nuclear and mitochondrial trees and between concatenated and coalescent-based inferences. Our results provide an alternative phylogenetic hypothesis to the existing mitochondrial DNA estimate of phylogeny. This new hypothesis provides a more accurate framework in which to explore trait evolution and examine the evolution of the mitochondrial genome in this group.

Spatio-temporal changes in the genetic structure of the Passerina bunting hybrid zone.

Although theoretical models predict that the structure of a hybrid zone can change under a variety of scenarios, only a few empirical studies of hybrid zones have unequivocally demonstrated zone movement. These studies are rare because few data sets exist that include repeated, temporally spaced, samples of the same hybrid zone. We analysed mitochondrial DNA haplotype data from samples separated by 40-45 years from across the Passerina amoena (Lazuli Bunting) and Passerina cyanea (Indigo Bunting) hybrid zone to investigate whether the genetic structure of this zone has changed during that interval. Both cline and generalized linear mixed modelling analyses uncovered a significant narrowing and a substantial westward shift of the Passerina bunting hybrid zone, clearly illustrating hybrid zone movement. The cause of the change may be due to a combination of ecological, demographic and behavioural factors. Our results predict that the width of the hybrid zone will continue to narrow over time, a finding consistent with reinforcement theory.

Reconstructing the phylogeny of "Buarremon" brush-finches and near relatives (Aves, Emberizidae) from individual gene trees.

Gene trees are often assumed to be equivalent to species trees, but processes such as incomplete lineage sorting can generate incongruence among gene topologies and analyzing multilocus data in concatenated matrices can be prone to systematic errors. Accordingly, a variety of new methods have been developed to estimate species trees using multilocus data sets. Here, we apply some of these methods to reconstruct the phylogeny of Buarremon and near relatives, a group in which phylogenetic analyses of mitochondrial DNA sequences produced results that were inconsistent with relationships implied by a taxonomy based on variation in external phenotype. Gene genealogies obtained for seven loci (one mitochondrial, six nuclear) were varied, with some supporting and some rejecting the monophyly of Buarremon. Overall, our species-tree analyses tended to support a monophyletic Buarremon, but due to lack of congruence between methodologies, resolution of the phylogeny of this group remains uncertain. More generally, our study indicates that the number of individuals sampled can have an important effect on phylogenetic reconstruction, that the use of seven markers does not guarantee obtaining a strongly-supported species tree, and that methods for species-tree reconstruction can produce different results using the same data; these are important considerations for researchers using these new phylogenetic approaches in other systems.

Historical divergence and gene flow: coalescent analyses of mitochondrial, autosomal and sex-linked loci in Passerina buntings.

Quantifying the role of gene flow during the divergence of closely related species is crucial to understanding the process of speciation. We collected DNA sequence data from 20 loci (one mitochondrial, 13 autosomal, and six sex-linked) for population samples of Lazuli Buntings (Passerina amoena) and Indigo Buntings (Passerina cyanea) (Aves: Cardinalidae) to test explicitly between a strict allopatric speciation model and a model in which divergence occurred despite postdivergence gene flow. Likelihood ratio tests of coalescent-based population genetic parameter estimates indicated a strong signal of postdivergence gene flow and a strict allopatric speciation model was rejected. Analyses of partitioned datasets (mitochondrial, autosomal, and sex-linked) suggest the overall gene flow patterns are driven primarily by autosomal gene flow, as there is no evidence of mitochondrial gene flow and we were unable to reject an allopatric speciation model for the sex-linked data. This pattern is consistent with either a parapatric divergence model or repeated periods of allopatry with gene flow occurring via secondary contact. These results are consistent with the low fitness of female avian hybrids under Haldane's rule and demonstrate that sex-linked loci likely are important in the initial generation of reproductive isolation, not just its maintenance.

Speciation in Passerina buntings: introgression patterns of sex-linked loci identify a candidate gene region for reproductive isolation.

Sex-chromosomes are thought to play an important role in speciation, but few studies of non-model organisms have investigated the relative influence of multiple sex-linked markers on reproductive isolation. We collected 222 individuals along a geographical transect spanning the hybrid zone between Passerina amoena and P. cyanea (Aves: Cardinalidae). Using maximum-likelihood cline fitting methods, we estimated locus-specific introgression rates for 10 z-linked markers. Although the cline width estimates ranged from 2.8 to 584 km, eight of 10 loci had cline widths between 224 and 271 km. We also used coalescent-based estimates of locus-specific divergence times between P. amoena and P. cyanea to test a recently proposed hypothesis of an inverse relationship between divergence time and cline width but did not find a significant association. The narrow width (2.8 km) of the cline estimated from the VLDLR9 locus indicates strong selection retarding introgression of alleles at this locus across the hybrid zone. Interestingly, a mutation in the very low density lipoprotein receptor (VLDLR) gene, in which VLDLR9 is an intron, is known to reduce the egg-laying ability of some chickens, suggesting a possible link between this gene region and reproductive isolation between P. amoena and P. cyanea. These results underscore the importance of sampling multiple loci to investigate introgression patterns across a chromosome or genome and support previous findings of the importance of sex-linked genes in speciation.

Haldane's rule in an avian system: using cline theory and divergence population genetics to test for differential introgression of mitochondrial, autosomal, and sex-linked loci across the Passerina bunting hybrid zone.

Using cline fitting and divergence population genetics, we tested a prediction of Haldane's rule: autosomal alleles should introgress more than z-linked alleles or mitochondrial haplotypes across the Passerina amoena/Passerina cyanea (Aves: Cardinalidae) hybrid zone. We screened 222 individuals collected along a transect in the Great Plains of North America that spans the contact zone for mitochondrial (two genes), autosomal (four loci) and z-linked (two loci) markers. Maximum-likelihood cline widths estimated from the mitochondrial (223 km) and z-linked (309 km) datasets were significantly narrower on average than the autosomal cline widths (466 km). We also found that mean coalescent-based estimates of introgression were larger for the autosomal loci (0.63 genes/generation, scaled to the mutation rate mu) than for both the mitochondrial (0.27) and z-linked loci (0.59). These patterns are consistent with Haldane's rule, but the among-locus variation also suggests many independently segregating loci are required to investigate introgression patterns across the genome. These results provide the first comprehensive comparison of mitochondrial, sex-linked, and autosomal loci across an avian hybrid zone and add to the body of evidence suggesting that sex chromosomes play an important role in the formation and maintenance of reproductive isolation between closely related species.

Phylogenetic conservatism and antiquity of a tropical specialization: army-ant-following in the typical antbirds (Thamnophilidae).

One of the most novel foraging strategies in Neotropical birds is army-ant-following, in which birds prey upon arthropods and small vertebrates flushed from the forest floor by swarm raids of the army-ant Eciton burchellii. This specialization is most developed in the typical antbirds (Thamnophilidae) which are divisible into three specialization categories: (1) those that forage at swarms opportunistically as army-ants move through their territories (occasional followers), (2) those that follow swarms beyond their territories but also forage independently of swarms (regular followers), and (3) those that appear incapable of foraging independently of swarms (obligate followers). Although army-ant-following is one of the great spectacles of tropical forests, basic questions about its evolution remain unaddressed. Using a strongly resolved molecular phylogeny of the typical antbirds, we found that army-ant-following is phylogenetically conserved, with regular following having evolved only three times, and that the most likely evolutionary progression was from least (occasional) to more (regular) to most (obligate) specialized, with no reversals from the obligate state. Despite the dependence of the specialists on a single ant species, molecular dating indicates that army-ant-following has persisted in antbirds since the late Miocene. These results provide the first characterization of army-ant-following as an ancient and phylogenetically conserved specialization.

Gene sampling strategies for multi-locus population estimates of genetic diversity (theta).

BACKGROUND: Theoretical work suggests that data from multiple nuclear loci provide better estimates of population genetic parameters than do single loci, but just how many loci are needed and how much sequence is required from each has been little explored. METHODOLOGY/PRINCIPLE FINDINGS: To investigate how much data is required to estimate the population genetic parameter theta (4N(e)mu) accurately under ideal circumstances, we simulated datasets of DNA sequences under three values of theta per site (0.1, 0.01, 0.001), varying in both the total number of base pairs sequenced per individual and the number of equal-length loci. From these datasets we estimated theta using the maximum likelihood coalescent framework implemented in the computer program Migrate. Our results corroborated the theoretical expectation that increasing the number of loci impacted the accuracy of the estimate more than increasing the sequence length at single loci. However, when the value of theta was low (0.001), the per-locus sequence length was also important for estimating theta accurately, something that has not been emphasized in previous work. CONCLUSIONS/SIGNIFICANCE: Accurate estimation of theta required data from at least 25 independently evolving loci. Beyond this, there was little added benefit in terms of decreasing the squared coefficient of variation of the coalescent estimates relative to the extra effort required to sample more loci.

Deletion of the fcgamma receptor IIb in thymic stromal lymphopoietin transgenic mice aggravates membranoproliferative glomerulonephritis.

Engagement of immunoglobulin-binding receptors (FcgammaR) on leukocytes and other cell types is one means by which immunoglobulins and immune complexes activate effector cells. One of these FcgammaRs, FcgammaRIIb, is thought to contribute to protection from autoimmune disease by down-regulation of B-cell responsiveness and myeloid cell activation. We assessed the role of FcgammaRIIb in a mouse model of cryoglobulin-associated membranoproliferative glomerulonephritis induced by overexpression of thymic stromal lymphopoietin (TSLP). TSLP transgenic mice were crossbred with animals deficient for FcgammaRIIb on the same genetic background (C57BL/6). Renal pathology was assessed in female and male animals (wild-type, FcgammaRIIb-/-, TSLP transgenic, and combined TSLP transgenic/FcgammaRIIb-/- mice) after 50 and 120 days, respectively. FcgammaRIIb-/- mice had no significant renal pathology, whereas overexpression of TSLP induced a membranoproliferative glomerulonephritis, as previously established. TSLP transgenic FcgammaRIIb-/- mice appeared sick with increased mortality. Kidney function was significantly impaired in male mice corresponding to aggravated glomerular pathology with increases in glomerular matrix and cellularity. This resulted from both a large influx of infiltrating macrophages and increased cellular proliferation. These results emphasize the important role of FcgammaRIIb in regulating immune responses and suggest that modulation of Fcgamma receptor activation or expression may be a useful therapeutic approach for treating glomerular diseases.